Phase Separation of Multi-Component Bose-Einstein Condensates of Trapped Atoms 
and Molecules with a Homonuclear Feshbach Resonance 
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We investigate phase separation of Bose-Einstein condensates (BECs) of two-component atoms 
and one-component molecules with a homonuclear Feshbach resonance. We develop a full model 
for dilute atomic and molecular gases including correlation of the Feshbach resonance and all kinds 
of interparticle interactions, and numerically calculate order parameters of the BECs in spherical 
harmonic oscillator traps at zero temperature with the Bogoliubov's classical field approximation. 
As a result, we find out that the Feshbach resonance can induce two types of phase separation. 
The actual phase structures and density profiles of the trapped gases are predicted in the whole 
parameter region, from the atom dominant regime to the molecule dominant regime. We focus on 
the role of the molecules in the phase separation. Especially in the atom dominant regime, the role 
of the molecules is described through effective interactions derived from our model. Furthermore we 
show that a perturbative and semi-classical limit of our model reproduces the conventional atomic 
BEC (single-channel) model. 

PACS numbers: 03.75.Mn, 05.30. Jp, 67.85.Fg 



I. INTRODUCTION 

Since the realization of Bose-Einstein condensation in 
alkali atoms in 1995 pHS, studies on cold atomic gases 
have greatly advanced Q. The cold atoms are pre- 
cisely controlled by electromagnetism and optics, and 
flexible to design quantum systems. In particular, inter- 
atomic interactions are effectively controlled in the Fes- 
hbach resonance technique [J], which realizes, e.g., artifi- 
cial collapse of atomic BECs @, BEC-BCS crossover in 
Fermi gases @ , and phase separation of multi-component 
BECs 0, [| related to this work. 

In this paper, we consider BECs of multi-component 
atoms. The components of the atoms are determined 
by their atomic species, mass numbers, and spin states 
and manually designed in actual experiments. In fact, 
two-component BECs of 87 Rb atoms with different spin 
states and three-component BECs of 23 Na atoms with 
spin triplet (spinor BECs) are produced in the end of 
1990s Physics of the multi-component BECs is 

one of the current issues in quantum physics. 

The multi-component BECs can exhibit immiscible 
phase separation, where the BECs do not exist as a mix- 
ture in equilibrium but produce domain structures made 
of different components of the BECs. Although such a 
phase separation is commonly seen in various liquids, the 
phase separation of dilute gases is a special feature of 
quantum systems dominated by energy in stead of en- 
tropy inducing the miscibility of classical gases. In fact, 
the phase separation of two-component atomic BECs 
is due to competition of the repulsive energies between 
same and different components of the BECs (Ea - llij . The 
phase separation reflects quantum statistics and interac- 
tions and is thus important in quantum physics. 

The behaviors of the atomic BECs, e.g., the phase sep- 



aration, can be controlled in the Feshbach resonance tech- 
nique. In this technique, molecules made of the atoms 
are introduced through the Feshbach resonance, which 
is a resonance between an interatomic scattering state 
(open-channel) and a molecular bound state with a dif- 
ferent spin structure (closed-channel). As a result, the 
molecules work as a resonant intermediary and effectively 
modifies the interaction in the open-channel. Energy de- 
tuning of the resonance can be controlled with external 
magnetic field owing to the Zeeman effect [15[. Recently 
two experimental groups have succeeded in inducing the 
phase separation in this technique: JILA group demon- 
strates it in BECs of 85 Rb and 87 Rb atoms with homonu- 
clear 85 Rb molecules [7|; Gakushuin group demonstrates 
it in the spinor BECs of 87 Rb atoms @ . 

The purpose of this paper is to describe the phase sep- 
aration and domain structures of the multi-component 
BECs of the atoms and molecules with the homonuclear 
Feshbach resonance. Here the molecules can play an ac- 
tive role and we are particularly interested in the rela- 
tionship between the molecules and phase structures. 

In order to accomplish the purpose, we develop a full 
model for the atoms and molecules including correlation 
of the Feshbach resonance, all kinds of interparticle inter- 
actions, and trap potentials. This model enables to pre- 
dict the actual phase structures and density profiles of the 
trapped BECs in the whole parameter region, from the 
atom dominant regime to the molecule dominant regime, 
and to investigate the overall behaviors of the molecules. 

In general, molecular formations in quantum gases 
can induce various phase structures. In quantum many- 
particle physics, molecules appear as a kind of many- 
body correlations, e.g., the molecule-pair (BEC-BCS) 
crossover for the molecule made of a pair of fermions [18| . 
The crossover theory is originally developed in the 
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fermionic systems and adapted in the boson-fermion en- 
vironments four-body correlations [2(j, and so on. 
In particular, one of the authors and his collaborators 
develop quasi-chemical equilibrium theory for diatomic 
molecules to the all kinds of quantum statistical envi- 
ronments including the bosonic environment related to 
the present work [2l| . In addition, the molecular forma- 
tions in the bosonic systems are realized in actual exper- 
iments [HI HI[ and studied in other theories [24| [H[ . 

In the atom dominant regime, the contribution of the 
molecules can be described as effective interactions for a 
partial atomic space in the full model, where the residual 
molecular space induces the effective interactions owing 
to the atom-molecule interactions and correlation. In 
general, the partial system and effective interactions are 
artificially defined from the full system and not unique. 
The theory of the effective interactions is a powerful tool 
to study physics [26| . 

A conventional approach to treat the Feshbach res- 
onance in the atom dominant regime is introduction 
of an effective interaction to the atomic BEC (single- 
channel) model in a perturbative way 0, EH, fl7| , where 
the molecules implicitly contribute to the atoms only 
through the effective interaction. In Sec IIV1 we show 
that a perturbative and semi-classical limit of our model 
reproduces the single-channel model. 

This paper is organized as follows. In section [ill we 
propose our model and derive formulas with some ap- 
proximations. In section IIII( we show a macroscopic 
analysis of the phase separation and numerical results of 
the phase structures and density profiles of the trapped 
gases. In section ITVl we define the effective interactions 
in our model to see the role of the molecules in the atom 
dominant regime. In addition, we discuss the correspon- 
dence relation and difference between our model and the 
single-channel model. Section [V] is devoted to summary 
and perspective. 



II. FORMULATION 

We consider the quantum gases of two-component 
bosonic atoms, labeled a and 6, and one-component 
molecules, labeled m, formed from pairs of the a atoms. 
Assume that the gases are trapped in vacuum, and total 
numbers, N at and Nm, of the a and b atoms are fixed and 
related to respective particle numbers, N a , N^, and N m , 
of the a atoms, b atoms, and molecules as 

N at = N a + 2N m , N bt = N b , (1) 

where N a and N m can vary through the Feshbach reso- 
nance. 

In this paper, the subscripts for the labels a, b, and m 
indicate the corresponding particles. In addition, we take 
H = 1 for the reduced Planck constant as a convention. 



The Hamiltonian is denoted by 



H = 


E 
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E E 


+ Hfk 




ex— a,fc,m 




a—a,b,m 0—a,b,m 





with the free-motion parts, H ai interparticle interaction 
parts, H a i3, and Feshbach resonance part, Hfr- 
The free-motion part H a in eq. ((2]) is defined as 



H a = J dr *l(r)h a {r)V a (r) (3) 

with the boson field operator *S> a (r) obeying the commu- 
tation relations, 

[* Q (r),^(r')] =0 (4) 

and 

[*a(r),*J(r')] =5 aP 5{r-r'). (5) 

The single-particle Hamiltonian h a (r) in eq. ([3]) is given 

as 

V 2 

K (r) = - — + V a (r) + M a (6) 

with the particle mass m a , trap potential V a (r), and 
mass term M a . In the nonrclativistic model, the values 
of the mass terms have no meaning except for the relative 
value between those of the atom and molecule. Thus we 
can put M a = M b = and M m = AE(r) + 2V a {r) - 
V m {r) with the energy detuning 

AE(r) = V m (r) + M m - 2(V a (r) + M a ) (7) 

without loss of generality. The energy detuning AE(r) 
can be varied by controlling the external magnetic field 
owing to the Zeeman effect. 

The interaction part H a p in eq. © is defined as 

H afi = 9 -f J dr ^(r)*+(r)^(r)* a (r) (8) 

with the contact-type interactions in the pseudo- 
potential approach. This approach must be valid for the 
low-energy scattering dominated by the s-wave scatter- 
ing length and also for the interactions in the cold dilute 
gases. For the BECs, we can put 

m a + mp 

9a/3 = 3fia = 27T a a/3 (9) 

m a mp 

with the bare s-wave scattering length a a p between the a 
and /3 particles in vacuum according to the regularization 
scheme. This approach has been widely applied to the 
cold atoms and molecules Q. Note that the interaction 
part in eq. (|8]) does not include effect of the Feshbach 
resonance correlation. 

We here supplement explanation of the regularization 
procedure for the contact-type interactions. In general, 
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the coupling constants, g a p, are theoretical dummy pa- 
rameters to reproduce relations among physical values 
through the regularization procedure. In fact, direct per- 
turbative calculation for g a p can produce the divergence 
difficulty owing to insufficient treatment of singularity of 
the two-body wave-functions. In the above paragraph, 
to deal with the difficulty, we introduce the renormalized 
coupling constants in eq. ^ with the two-body wave- 
functions in vacuum. That is because the singularity of 
the two-body wave-functions is assumed as a extremely- 
short-range property separated from the global property 
described in the asymptotic fields for the BECs. 

The Feshbach resonance part -Hfr in eq. ([2]) is defined 

as 



Hfr = A / dr 



*LW(^a(r)) 2 + (*t(r)) 2 vI/ m (r)' 



( 10 ) 

also in the pseudo-potential approach, where we apply 
the same regularization procedure for the BECs. The 
absolute value of the A in eq. (|10p can be obtained from 
the detuning dependence of the s-wave scattering length 
with the Feshbach resonance in vacuum, 

fill 

ia 2tt AE' [ ' 

as shown in appendix^] In addition, we can freely deter- 
mine the sign of the A in eq. (fTU]l because it is related to 
the phases of the field operators. The a aa (AE) and a aa 
in eq. (|11[) can be obtained in scattering experiments in 
principle. No free parameter thus survives in this model. 

Here we apply the Bogoliubov's classical field approx- 
imation [29[ . First we introduce BECs defined with the 
BEC order parameters, <j) a (r) = (* Q (r)), by using the 
Bogoliubov's method, ^ a (r) — (f) a (r) + & a (r). Next we 
neglect the thermal and correlational part, 



*a(r) « a (r), 



(12) 



in the zero-temperature equilibrium. This approxima- 
tion must be valid in BEC dominant regime as shown 
in the actual experiments 0, HI , and can not be applied 
to strongly-correlated Bose gases with little or no BECs. 
The latter situation is beyond the purpose of this paper. 

The order parameters, <J) a (r), in the zero-temperature 
equilibrium are obtained from minimization of the energy 
given in eq. @ with eq. (|12[) under the number constraint 
m eq. ((TJ with 



dr\<f> a (r)\ 



(13) 



In the Lagrange's multiplier method, we obtain the cou- 
pled non-linear Schrodinger-type equations, 

lla,4>a = ( K(r) + ^ ffaa|</>a| 2 Wa + 2\(j) m (j)* a , (14) 
a.=a,b.m 



fi b <j) b = h b (r) + 



E 

a— a,6,m 



9bo 



(15) 



and 



/in 



= h m (r) + 



E 

a—a,b,m 



Quia tpot J + A#, (16) 



and chemical equilibrium equation, 

— f-^rn i 



(17) 



for the chemical potentials, \i a . 

Lastly we summarize the treatment of BEC in the 
present approach. In principle, BEC is a phase transition 
defined in macroscopic systems; On another front, the 
confinement traps produce microscopic systems, where 
the phase transition is not well-defined owing to the non- 
negligible finite size effect. In order to introduce BEC 
to the microscopic system, we must determine the treat- 
ment of the finite size effect. In the present approach, we 
firstly apply the Bogoliubov's method in eq. (fT2j) to in- 
troduce BEC without the finite size effect, and afterward 
introduce the traps and finite size effect. As a result, 
the BECs are defined in the microscopic system with the 
finite size effect. In fact, the phase separation of the 
BECs appears not as a sharp density separation but as 
a smooth density crossover as demonstrated in the next 
section 1301 • 



III. PHASE SEPARATION AND DOMAIN 
STRUCTURES 

In this section, we show results of the calculation ac- 
cording to eqs. (|13l) - (fl7|) . Parameters are defined in sub- 
section IIII Al Calculational scheme is given in subsec- 
tion IIIIBI We show general property and typical behav- 
iors of the density profiles in subsection IIII CI and phase 
structures in subsection IIII Dl 



A. Parameters 

In order to demonstrate the calculation, we assume a 
situation as follows: First the a and b atoms have mostly 
the same masses, m a = mi, = m, and coupling constants, 



9a 



g a b = gbb = g. 



(18) 



Second the atoms and molecules are trapped in the spher- 
ical harmonic oscillator potentials denoted by 



V a (r) = )^m a w 2 r 2 



(19) 



with the trap frequency u>. Here m m &s 2m for the nonrel- 
ativistic diatomic molecule, and the spatial dependence 
of the energy detuning AE(r) in eq. (0 is negligible. 
This assumption corresponds to the case of atoms with 
two spin-components, and leads to a typical demonstra- 
tion of the phase separation. 
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We also assume the hard-sphere hypothesis: The 
atoms and molecule have hard-sphere with radii R a 
(= Rb) and R m , respectively, and the scattering lengths 
in eq. © become 

a a p = R a + Rp- (20) 
According to eq. ©, we then obtain 

3 

9a,n = 9bm = g (1 + «) 9 (21) 

and 

9mm = (22) 

with k = Rraj R a , where g a p can be described only with 
the two positive parameters g and k. Note that this 
hypothesis is just for convenience to fix the parameters 
and simplify the model for demonstration. In principle, 
the scattering lengths should be determined in scattering 
experiments. 

We summarize the parameters. In the harmonic oscil- 
lator units, m = 1, lj = 1, and H = 1 without loss of 
generality, and the remaining parameters are g, k, A, Nb, 
N at (— N a + 2N m ), and AE. Note that the number ratio 
of the molecules, N m /N at , is a monotonically decreasing 
function of the detuning energy, AE, and we can then 
take N m as a parameter instead of AE. 

B. Calculational scheme 

In order to compute eqs. (fl^|) - (jT5|l . we utilize the 
imaginary-time relaxation method with the spatial dif- 
ference scheme. In this method, we successively improve 
a trial state starting from an input initial state to obtain 
a solution. Thus the resulting solution may generally de- 
pend on the initial state if several solutions coexist in a 
numerical precision. 

In this paper, we give the initial state in the follow- 
ing scheme. First we solve them with an arbitrary ini- 
tial state at N m ~ 0, where we can not find the initial 
state dependence. Next we successively solve them with 
asymptotical increment of N m by substituting the ob- 
tained state with tiny numerical noises to avoid numer- 
ical instability for the next initial state. We iterate this 
process until N m reaches a given number. 

In addition, we also perform the same scheme but 
started from N m ~ N at /2 with asymptotical decrement 
of N m . As a result, we confirm that both results agree 
in the numerical precision. 

C. General property 

The phase separation of the BECs is due to a discon- 
tinuous phase transition in terms of a ratio between two 
of the component densities. This phase transition can be 



described as instability for local density fluctuations in 
the uniform mixture. In general, the instability indicates 
not only the phase separation - appearance of lower lim- 
its of the densities in equilibrium - but also collapse - 
appearance of the divergent densities. If a system does 
not induce the collapse, the instability is directly linked 
to the phase separation. 

In the present system, the instability indicates the 
phase separation because this system has no attraction 
to induce the collapse. In fact, the only attraction in this 
system is the Feshbach resonance part in eq. (|10[) . which 
can not induce the collapse because the power of the den- 
sity dependence of the attraction energy in eq. (|TU)) . 3/2, 
is less than that of the repulsion energy in eq. (JSJ, 2. 

The stability condition at the zero temperature can be 
described as minimum condition of the energy density 
functional with the number constraints. The detail of 
the stability condition is summarized in appendix |B| 

The stability condition suggests that the three- 
component mixture is always unstable except for A = 
and breaks into immiscible domains. In other words, this 
system exhibits the phase separation whenever the Fes- 
hbach resonance correlation in eq. (fTO"]) exists. This re- 
sult depends on the assumption for the atoms and hard- 
sphere hypothesis in eqs. (|2~T|) and (f22j). and is indepen- 
dent of values of the parameters in subsection IIII Al 

Moreover, according to the stability condition, possi- 
ble phases in the immiscible domains are limited to (6), 
(m), and (a, m), where the labels in the round brackets 
represent the components of the existing particles in the 
equilibrium. The phase (b) is always stable or metastablc 
while the other phases have their own stability conditions 
described in appendixlB] In particular, the phase (m) has 
a lower limit for the molecular density n m and is stable 
or metastable for the larger n m than the limit. 

In the microscopic trap system, the phase separation 
does not make the sharp domain structure except for 
the theoretical Thomas-Fermi limit, N a t,Nb —> oo, be- 
cause of the finite size effect. In the actual experiments, 
the numbers of the atoms are limited below 10 4 ~ 6 in 
the evaporative cooling method, and ranges of interfaces 
between the domains, i.e., the healing lengths, are not 
negligible. 

Furthermore the above macroscopic analysis for the 
phase separation does not give any information for pro- 
files and compositions of the domains. That is because 
they depend on the interface effects, which can not be 
given as extensive thermodynamical variables. In order 
to obtain the actual domain structure, the microscopic 
analysis is needed beyond the macroscopic analysis. 

In fig. [TJ we show typical behaviors of the density 
profiles in the spherical system with r = \r\ when 
g = 0.1 x mT 3 ! 2 ^- 1 ! 2 , k = 2, A = 0.1 x m^u; 1 / 4 , 
N at = N b = 5x 10 4 , and N m = 5 x 10 1 (a), 5 x 10 2 (b), 
and 5 x 10 3 (c). If we change the values of the parameters 
as in the next subsection, the qualitative behaviors are 
mostly similar to those in fig. [TJ 

Fig. [TJl represents the profiles in the atom dominant 
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FIG. 1: (Color online) The density profiles when 



10 



0.1 x 



-1/2 



, k — 2, A = 0.1 xm 



"3/4,. 



1/4 ,iV at = iV, ) = 5xl0 4 



and iV m = 5 x 10 1 (a), 5 x 10 2 (b), and 5 x 10 3 (c). The solid, 
dashed, and dotted lines denote those of the a atoms, b atoms, 
and molecules, respectively. 



regime, N rn /N at = 1(T 3 . Although N m /N at a 0, the 
molecules affect the atomic distributions, which are ob- 

l 1 2 

viously different from the symmetric ones, </> a (r) — 

l 1 2 

</>f,(r) , at the N m = limit. It is due to the Fesh- 
bach resonance part in eq. ([TO)) and suggests the phase 
separation. 

The phase separation is more clearly seen in fig. [TJd 
for N m /N at = 10~ 2 . As N m increases, the b atoms are 
gradually localized outside the a atoms and molecules. 
According to the macroscopic analysis, these profiles in- 
dicate an unsharp double-domain structure made of the 
two phases (a, m) and (b) in common with those in fig.[TJi. 

Fig. [Ip represents the density distributions in the 
atom-molecule coexistence regime, N m /N at = 10 _l . As 
N m increases, the a atoms are gradually localized in the 
middle region between the inner molecules and outer b 
atoms. According to the macroscopic analysis, these pro- 
files indicate an unsharp triple-domain structure made of 
the three phases (m), (a, m), and (b). Here we empha- 
size that this structure can not be seen in the conven- 
tional single-channel model, which does not include the 
molecular phase (m) and must be valid only in the atom 
dominant regime. 



The constitution of the unsharp domain structures is 
due to minimization of the energy made up of the kinetic 
and trap parts in eq. ((3|), interaction parts in eq. ((8|), 
and Feshbach resonance part in eq. (JTUJ) . The domain 
profiles principally reflect the surface and interface ener- 
gies derived from the kinetic parts. The domain arrange- 
ment principally reflects the energy competition among 
the trap, interaction, and Feshbach resonance parts. In 
a rough explanation, the a atoms tend to take the cen- 
ter position in the atom dominant regime (figs.QJ, and b) 
because of the attraction in the Feshbach resonance part, 
which gives an advantage to enlarge the densities of the 
a atoms and molecules; the molecules tend to take the 
center position in the atom-molecule coexistence regime 
(fig. QJ) because of the mass difference in the trap poten- 
tials in eq. (TT^|) and oscillator lengths, 1/ ^m a u. 



D. Phase structures 

As demonstrated above, the phase separation induces 
the unsharp domain structures in the microscopic system. 
Some of the density profiles can take a hollow structure 
in the spherical system as shown in fig. [TJ 

In order to describe phase structures of the density 
profiles, we introduce indicator quantities defined as 





0a (0) 


2 


max 


|<M 


\l 2 " 

r)| 



(23) 



for a — a, b, and m. For the non-hollow-type profile, 
the maximal density is at the center, r — 0, and v a = 1; 
otherwise v a < 1. 

In fig. [21 we plot the phase structures for the a atoms 
(a) and b atoms (b) with the contour lines for the v a 
in eq. (|23[) . where we omit to plot for v m because the 
molecules do not exhibit the hollow-type profile and 
v m = 1. Here we choose the same values in fig. [1] for 
g, k, and A and give the total number of the atoms as 
N a t + Nb = 10 5 . The cases in fig. [1] correspond to those 
on the N a t — A& = line in fig. [21 Note that edges of 
the variable area in fig. [5] are singular, where one of the 
particle numbers is zero. 

Fig. [5] represents the smoothly-varying phase struc- 
tures. In the smaller N m regions than the v a = 0.99 
lines, v a w 1 and the corresponding particles take the 
non-hollow-type profile. In the larger N m regions than 
the v a = 10~ 2 lines, v a w and the corresponding parti- 
cles completely take the hollow-type profile. The value of 
v a is always larger than that of v . It indicates that, as 
N m increases, the double-domain structure appears first 
and afterward crosses over into the triple-domain struc- 
ture. As Nb increases, the effect of the molecules and a 
atoms on the b atoms relatively decreases, and the densi- 
ties of the molecules and a atoms also decrease; Then the 
appearance of the domain structures delays as a result. 

Lastly we comment on the dependence on the param- 
eters fixed in fig. O First, as N at + Nb and g increase, 
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Nm / Nat 



FIG. 2: (Color online) The phase structures of the density 
profiles of the a atoms (a) and 6 atoms (b) when g = 0.1 X 
m -3/2 w -i/2 ) K = 2) > = o.l x m- 3/4 o; 1/4 , and 7V at + N b = 
1.0 x 10 5 . The solid, dashed, dot-dashed, and dotted lines 
denote the contour lines for v a = 0.01, 0.1, 0.5, and 0.99, 
respectively. 

the system approaches to the Thomas- Fermi limit, where 
the domain structures appear sharply. Second, as |A| 
increases, the a atoms and molecules tend to mix and 
concentrate owing to the attraction in the Feshbach res- 
onance part, and the domain structures can easily appear 
owing to the concentration of the a atoms and molecules 
according to stability conditions in eqs. (|B7[) and (IB8|) . 
Third k does not affect the qualitative properties of the 
phase structures as confirmed in the stability conditions 
and also in the numerical analysis for various values of k. 

IV. EFFECTIVE INTERACTIONS 

In this section, we discuss the role of the molecules in 
the phase separation and double-domain structures in the 
atom dominant regime by introducing the effective inter- 
actions. In addition, we show that a perturbative and 



semi-classical limit of the present approach reproduces 
the single-channel model. Validity of the single-channel 
model is also discussed with numerical demonstration. 

In subsection IIVA1 we define the effective interactions 
for the present model. In subsection IIVB1 we take the 
limit and obtain the single- channel model. In subsec- 
tion UvUl we discuss the role of the molecules. In sub- 
section IIVD1 we compare the present and single-channel 
models. 



A. Exact formulation 

In order to define the effective interactions, we focus 
on a partial system only consisting of the BECs of the a 
and b atoms, where the number of the a atoms, N a , in 
the partial system is counted without the atoms in the 
molecules and does not agree with the total number in 
the full system, N at ^ N a . 

Since the partial system does not include the 
molecules, the BEC order parameters should obey 

fra{r)<t>a = (h a (r)+ ^2 9aa{r)\<t) a ^Ha (24) 

^ a=a,b ' 

and 

fo{r)4>b = (h b {r)+ ^2 9ba(r)\(f> a \ 2 )<fib (25) 

^ a=a,b ' 

in the same manner in section [TTl where the effective 
chemical potentials, /2 a (r), and effective coupling con- 
stants, g a p{r), can have the spatial dependence in gen- 
eral. 

Exact definitions of the effective quantities must not 
affect physical values, e.g., the order parameters and par- 
ticle numbers, because physics must not depend on the 
artificial definitions. 

In order to obtain the exact definitions, we utilize the 
solutions of the full model in self-consistent framework, 
where the effective interactions can reproduce the origi- 
nal solutions obtained from eqs. (fT3|) - ([T7|) . 

The effective chemical potentials, ji a (r), in eqs. 
and (|25[) are defined as 

l 1 2 

p, a (r) = n a - g am \(/> m (r)\ (26) 

with the order parameter 4> m (r) calculated from the full 
model, where fi a (r) reflects the background energy de- 
rived from the atom-molecule interactions and indicates 
renormalization of the mass-energy owing to the effective 
interactions. 

The effective coupling constants, g a p{r), in eqs. (|24|) 
and (|25[) are defined as <? a/ 3(r) = g a p except for 

_ 47ro aa r cj> m (r) <f>; n (r) 1 
9aai) = ~m^~ [(Mr)) 2 + mrW\ ^ 

with the order parameters 4>a( r ) and <t> m { r ) calculated 
from the full model, where g a ,a{i') is renormalized with 
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the Feshbach resonance correlation and produces one 
of the effective interactions. According to eq. (fTS]). we 
rewrite eq. ([27]) as 



laa(r) 



4na n 



2A 2 



(28) 



with 



Mm( r ) = Mn 



(h m {r)4> m {r)) 
4> m (r) 



- ^2g ma \<i)a{r)\ 2 ■ (29) 



Furthermore, according to cqs. (fTTj) and (|T4"|). the /i m in 
eq. (|29l) can be rewritten as 



Mm — 2 



(WW) 

Mr) 



(30) 



The exact definitions are applicable to the whole pa- 
rameter region and reproduce the phase structures in the 
previous section. In fact, the double-domain structure is 
principally derived from g aa (r); The triple-domain struc- 
ture reflects all of the effective quantities. 

The effective interactions provide a physical interpre- 
tation from the viewpoint of the partial system when 
we obtain the actual values of the effective quantities, 
which are determined from the full system, not the par- 
tial system, except for some trivial situations, e.g., the 
decoupling or perturbative limit. 



B. The single-channel model 

We here take the perturbative and semi-classical limit 
of the present approach in the atom dominant regime, 
N a t ~ N a and N m /N at ~ 0, where the perturbation is 
started from vacuum of the atoms and molecule. 

In the limit, we directly obtain N at ~ N a and /i a (T*) ~ 
fi a in eq. (f2l)]) because 4> m (r) ~ 0. We thus focus on the 
g aa (r) in eq. (H7J). 

We also obtain the following relations: In the pertur- 
bative limit and not so close the resonance, 



la/3 0/3 







(31) 



because 4'p{ r ) ~ near vacuum. In the semi-classical 
limit, 



(h a (r)4> a (r)) 
<f> a (r) 



V a (r) + M a 



(32) 



owing to eq. ([5]) and the zero momenta of BECs, where 
the zero-point kinetic energies are neglected. 

By substituting eqs. (|3"Tj) and (|32|) into eq. (l28l) with 
eqs. ©, ([Ml), and flSU]), we obtain 



zaa 



4na a 



2A 2 
AE' 



(33) 



If V m (r) ss 2V a (r) as assumed in eq. (|T9|) . the AE in 
eq. (|33p does not have the spatial dependence. 

The effective interaction g a s ^ in eq. (|33j) agrees with 
that in the single-channel model, where the resonance 
changes a aa into a aa (AE) in eq. (fTI]) and g aa into gffl 
in eq. (1331) according to eq. ©. 

Except for the limit, the results in the single-channel 
model disagree with those in the full model. We label 
the quantities in the single-channel model with the su- 
perscript (sc) to distinguish them from the exact ones. 

In general, the perturbative approach must be valid 
when the perturbation is small and the nonperturba- 
tive state is appropriately selected. The former condi- 
tion does not matter here because of the regularization in 
eq. ([33]) , which corresponds to partial sum of the ladder- 
type perturbation series; The latter condition must be 
verified by comparing the perturbative approach with the 
exact approach. 



C. The role of the molecules 

In the atom dominant regime, the phase separation 
occurs between the a and b atoms and thus makes the 
double-domain structure not only in the full system but 
also in the partial system. 

In the macroscopic analysis for the partial system, the 
stability condition suggests the phase separation if 



9aa {r)gbb ~ g a b9ba < 0. 



(34) 



It is obtained in the same manner for the full system and 
becomes 



9aa(r) < g 



(35) 



in the situations in the previous section. Eq. (I35|) is al- 
ways true and indicates the phase separation because the 
first term in the right-hand side in eq. (|28[) equals to q 
and the second term is negative in the equilibrium [31] . 
This result is consistent with that for the full system. 

The role of the molecules in the phase separation is 
thus given in g aa (r) in the atom dominant regime, where 
the minority molecules hardly affect the other effective 
quantities, i.e., fla(f) ~ fJ. a , and behave as the reso- 
nant intermediary. It is clearly described in eq. (I33[) in 
the single-channel model limit and generally described in 
eq. (|28[) including the medium effect in the equilibrium. 
The difference between the results in the single-channel 
and full models is due to the medium effect. 



D. Comparison 

We here show the difference in the same situation in 
the previous section: g = 0.1 x m _3//2 cj _1 / 2 , n = 2, A = 
0.1 x m- 3 / 4 ^ 1 / 4 , and N at = N b = 5 x 10 4 . 

In fig. [3J we plot the density profiles in the single- 
channel and full models for N m /N at = 10~ 2 . The profiles 



8 



3 




1/2 

r I (H/moi) 



FIG. 3: (Color online) The density profiles when g = 0.1 x 
m -3/2 w -i/2 ) K = 2 , A = 0.1xm- s/4 uj 1/4 , N at = N b = 5xl0 4 , 
and iV m = 5 x 10 2 . The solid and dashed lines denote those 
of the a and b atoms in the present model; the thin-solid and 
thin-dashed lines denote those of the a and b atoms in the 
single-channel model, respectively. 



0.95 
1 

o 

TO 

- — - f0 

0.9 





' (a) 


v. 


■ — , 




a atom 


■ — 




b atom 









k 




(b) 


\ " 






\ 

\ 

\ 

\ 

\ 

\ 






— e>/g 

§aa(0)/g 




\ 

\ 

i 







N 



0.05 
m/Nat 



0.1 



in the full model are same as those in fig. [Tp. Note that 
the AE in eq. ([55)) is determined as a function of N m /N at . 

Although the domain structures are similar in the two 
models, fig.[3]represents obvious difference in the profiles. 
The difference is mainly caused by the medium effect 
ignored in eq. ([51]) because this case is near the semi- 
classical limit and not so near the perturbative limit. In 
fact, if we introduce the effect ignored in eq. ([31]) to the 
single- channel model, the improved results almost agree 
with those in the full model. 

In order to evaluate fidelities of the order parameters 
in the single-channel model, we define 



F a =± 
N Q 

for a — a and 6, where 



Nn 



dr U a (r)\ 



dr r a (r)^ c \r) (36) 



dr Ui sc '(r)| 2 . (37) 



The F a in eq. ([3"6| approaches to one when <pa°\f) and 
4> a (f) are similar. 

In fig. 0^, we plot the F a and F b in eq. ([3"o) , which 
gradually decrease as N m increases. Here Ft, is always 
smaller than F a because the fidelities reflect the volume 
factor and the b atoms occupy the outer side. 

In fig. [4p, we plot the g aa (0) in eq. ([58) an d 9aa in 
eq. ([3"3"[) , where the difference between them gradually in- 
creases as N m increases in consistency with fig. 0^. Here 
g aa (0) is smaller than <?i s a c ^ because of the medium effect 
shown in fig. [3] In fact, if we assume eq. (f32[) . <j) m {r) ~ 0, 
and g a a{r) ~ <?, the fJ>' m (r) in eq. (1551) becomes 



~AE- 



3k - 13 



-9 



\(f>a(r)\ + \4>b(r)\ 



(38) 



and then —fx' m (r) < AE and 3aa(^) < 



_(sc) 
,9a a 



for 



FIG. 4: (Color online) The fidelities (a) and effective coupling 
constants (b) when g = 0.1 x m~ 3//2 uj~ 1/2 , k = 2, A = 0.1 x 
m -s/4 w i/4 ) and N at = N b = 5x 10 4 . The solid and dashed 
lines denote F a and Fb in fig. [4^,, respectively, and gi^Vfl 1 an d 
g aa (0)/g in fig. Hp. 



The spatial dependence of g a a(f) does not play any 
active roles here [Hj]. In fact, if we use g aa (0) for the 

single-channel model instead of g^ in fig. El the im- 
proved results mostly reproduce the original profiles and 
the corresponding fidelities take almost one (> 0.999). 



V. SUMMARY AND PERSPECTIVE 

In this paper, we have studied the phase separation of 
the atomic and molecular BECs with the homonuclear 
Feshbach resonance by developing the full model formu- 
lated in section [TTl 

The Feshbach resonance can induce the phase sepa- 
ration, which produces the unsharp domain structures 
in the trap system as demonstrated in section Mil the 
double-domain structure in the atom dominant regime 
and triple-domain structure in the atom-molecule coex- 
istence regime. Here the molecules behave not only as the 
resonant intermediary but also as another BEC different 
from the atomic BECs. 

In the atom dominant regime, the role of the molecules 
is limited to the resonant intermediary and revealed 
through the effective interactions as discussed in sec- 
tion IIV1 Furthermore the effective interactions provide 
the single-channel model in the perturbative and semi- 
classical limit. The validity of the single- channel model 
depends on the medium effect. 

We here comment on the recent experiments of the 
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artificial phase separation of BECs in the Feshbach res- 
onance technique 0, H| • The experiments are performed 
in the atom dominant regime, and the phase separation 
refers to that of the atomic BECs in analogy with the 
double-domain structure in the present work. The triple- 
domain structure has not been realized in actual exper- 
iments. We are expecting further development of the 
experimental study in this field. 

The actual experiments are usually performed in 
anisotropic traps, not in spherical traps as in the present 
work. The asymmetry of the traps can induce anisotropic 
domain structures as shown in the recent experiments 0, 
8]. In principle, the anisotropic case can be described 
in the same framework in this paper by adopting the 
anisotropic traps as V a (r) in cq. (J6)). 

We also comment on the previous study on the quasi- 
chemical equilibrium theory for the BECs [2lJ . The con- 
cept and details of the work are quite different from those 
of the present work: the macroscopic effective theory vs. 
the microscopic model, the heteronuclear molecules vs. 
the homonuclear molecules, the negligible correlation vs. 
the resonant correlation, and so on. However two signifi- 
cant analogies are seen in the results in both works. First 
the atom-molecule phase separation appears as the coex- 
istence phases in the previous work and triple-domain 
structure in the present work. It is due to the interpar- 
ticle interactions and reflects the repulsion between the 
atomic and molecular BECs. Second the interaction en- 
ergies can contribute to the results as the renormalized 
energy detuning in the previous work and medium effect 
in the present work. 

In the macroscopic analysis, some additional works 
have recently been published. In ref. [27) . the quasi- 
chemical equilibrium theory is developed to the case near 
the heteronuclear Feshbach resonance with the resonant 
correlation. In ref. 28j, the homonuclear molecule case 
is studied for the one-component atoms. 

We here emphasize that the macroscopic analysis can 
not give any information on the actual domain structures. 
In fact, e.g., the occurrence of the coexistence phases in 
the quasi-chemical equilibrium theory merely indicates 
the possible existence of the triple-domain structure and 
does not determine the appearance and profile, which are 
determined in the microscopic analysis as demonstrated 
in this paper. 

Finally we notice the application limit of the present 
approach in terms of the validity of the Bogoliubov's clas- 
sical field approximation, which must be valid in the BEC 
dominant regime and not be applied to the strongly cor- 
related Bose gases. If the resonant correlation is quite 
strong and the situation is near the molecule dominant 
regime, the atoms may behave as the strongly correlated 
gas in analogy with the nondegenerate fermionic atoms 
in the BEC-BCS crossover. This situation is beyond the 
purpose of this paper and should be studied in another 
paper. 
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Appendix A: Feshbach resonance 

Let us consider elastic scattering between two atoms 
in vacuum with a Feshbach resonance. The Schrodinger 
equation can be written as coupled channel equations, 

(E-PHP)\yr) = (PHQ)\V%), (Al) 

{E-QHQ)\*%) = {QHP)\* P E ), (A2) 

for = P\^> E ) and = Q\$>e) with the projec- 

tion operators P (= P' = P 2 ) for the open-channel and 
Q (= = Q 2 ) for the closed-channel obeying P + Q = 1 
and PQ = QP = 0. 

The open-channel Hamiltonian PHP = H Q in eq. (|A1|) 
is given as 

H Q = K Q + V (A3) 

with the kinetic part K a of the relative motion in the 
center-of-mass system and the bare interaction part V 
denoted by the contact-type interaction, 

K{ct>q\V \<l>p) =g, (A4) 

where we introduce the arbitrary normalization factor 0t 
and energy-momentum eigen-state, 

Ko\<t>p) = ^\4>p) = E\4> p ), (A5) 

with the reduced mass m. 

In the regularization in eq. (j9]), we renormalize 

ffH> 1^ * [™Pty p \ Vo \i, E ) = HZ^° (A6) 
e^o J An x ■ Pl 1 ' m 

with the solid angle ftp and two-body scattering wave- 
function, 

where the bare s-wave scattering length eto is related to 
g as 

^^(l + ffA)- 1 (A8) 
with the divergent integral 
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The closed-channel Hamiltonian QHQ = H c in eq. (|A2 

gives energy eigen-states in the Q-space, 



H c \xn) = £n\Xn), 



(A10) 



including the resonant molecular state \xm) and energy 
detuning, em — AE, for n = M. 

According to eqs. (|A1|) and (|A2I) . we obtain 



(All) 



(E-H )\^ P E )=V E \^ P E ) 
with the effective interaction 

1 



V E = (PHQ) 



where we assume 



E - H c +ie 



(QHP) , (A12) 



V E « (PHQ) l XM l^l\ (QHP) (A13) 
E — AE + te 

owing to the complete set in the Q-space, 

Y / kn)(Xn\=Q, (A14) 



According to eqs. (|A"8|) and (|A18I) - (|A20|) . we obtain 

i e 



a = a 



2ir AE - S 



and 



with the resonance shift 



2A 2 A 



and width 



1 - 47ra A 



e ~ AA' 



(A21) 



(A22) 



(A23) 



(A24) 



which are independent of AE. In the second-order per- 
turbation of A, we obtain m/A-E <C 1 and eq. (fTTj) , 



Appendix B: Stability conditions 



and the resonance, \E — £ n /M S> \E — AE\. 

The resonance part is denoted by the contact-type cor- 
relation, 

Vm(x M \QHP\(f>p) = V2X, (A15) 

where the coefficient \[2 is due to symmetrization for 
the identical atoms. Note that the exchange part, 
|</>_p)/v2, does not contribute to the molecule forma- 
tion in eq. (|A15|) . 

In the regularization in eq. (|11[) . we renormalize 



(A16) 



where the two-body wave-function in the open-channel is 
written as 



E-H Q -V E +ie 



(v + Ve) \4> p ), 

(A17) 

with the dissipation factor into the molecule accord- 
ing to eq. (|Alip . The A and A in eq. (|A16j) obey 



A = A(l+gA)- i 
with the effective coupling constant 



9 = 9 



2A 2 
AE' 



(A18) 



(A19) 



The s-wave scattering length ao in the open-channel is 
written as 



a = |-(l + .gA)- 1 



(A20) 



The energy density £[n a , n^,, n m ] for the uniform BECs, 
f H) = J dr £, is given as 



£ = 



E\ ^ 9a/3 



q — a , b , rn {3=a,b,m 



2\\\n a n]l 2 (Bl) 



with the g a p (= gp a ) in eqs. (0, ([2"T]l. and (l2"2"j) . 

The phase (a, 6, to) is unstable because the Hessian 



n (abm) _ 



<9 2 £ [n Q ,n6,n TOJ 



dn a dnp 



(B2) 



for a, /3 = a,b,m is negative: 



<0. 



(B3) 



The phase (6, m) is unstable because the Hessian 
r <9 2 £[n a = 0, n fc ,n„ : 



•H (bm) = det 



for a, /3 = 6, to is negative: 



dn a dnp 



H (bm) = 2 

64 y 



K 

9 



7A 32 

K 81 



< 0. 



(B4) 



(B5) 



The phases (a, b) and (a) are unstable because of the 
asymptotic behaviors near n rn — 0. 

The phase (a, to) is stable or meta-stable when the 
Hessian 



H {am) = dot 



d 2 £[n a ,n b = 0,n r , 



dn a dnp 



(B6) 
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for a, p = a, m is positive: 



2?7 



1/2 



n a 3(1 + k) 
7i*n 2 



The phase (6) is always stable or meta-stable. 



The phase (m) is stable or meta-stable when 



> 0. 
(B7) 



16 A 



n 2 



3 (l + n)g 

because of the asymptotic behavior near n a = ?i& = 



(B8) 
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The present approach is merely a theory to give an ef- 
fective model for the microscopic system. As another ap- 
proach, we can apply the variational method with a trial 
state denoted by the coherent state for the bosons. Then 
we obtain the same formulation after all. 
The effective scattering length in the equilibrium can not 
show the drastic sign reversal at the resonance, differently 
from that in vacuum in eq. (|ll|l . because the resonant 
correlation is always attractive in the equilibrium. 
Here g aa (r) is a monotonically increasing function for the 
radial coordinate r and approaches to g at large r. 



